Transposon dynamics in the emerging oilseed crop Thlaspi arvense.

Genome evolution is partly driven by the mobility of transposable elements (TEs) which often leads to deleterious effects, but their activity can also facilitate genetic novelty and catalyze local adaptation. We explored how the intraspecific diversity of TE polymorphisms might contribute to the broad geographic success and adaptive capacity of the emerging oil crop Thlaspi arvense (field pennycress). We classified the TE inventory based on a high-quality genome assembly, estimated the age of retrotransposon TE families and comprehensively assessed their mobilization potential. A survey of 280 accessions from 12 regions across the Northern hemisphere allowed us to quantify over 90,000 TE insertion polymorphisms (TIPs). Their distribution mirrored the genetic differentiation as measured by single nucleotide polymorphisms (SNPs). The number and types of mobile TE families vary substantially across populations, but there are also shared patterns common to all accessions. Ty3/Athila elements are the main drivers of TE diversity in T. arvense populations, while a single Ty1/Alesia lineage might be particularly important for transcriptome divergence. The number of retrotransposon TIPs is associated with variation at genes related to epigenetic regulation, including an apparent knockout mutation in BROMODOMAIN AND ATPase DOMAIN-CONTAINING PROTEIN 1 (BRAT1), while DNA transposons are associated with variation at the HSP19 heat shock protein gene. We propose that the high rate of mobilization activity can be harnessed for targeted gene expression diversification, which may ultimately present a toolbox for the potential use of transposition in breeding and domestication of T. arvense.

Characterization of the active site in the thiocyanate-forming protein from Thlaspi arvense (TaTFP) using EPR spectroscopy.

Glucosinolates are plant thioglucosides, which act as chemical defenses. Upon tissue damage, their myrosinase-catalyzed hydrolysis yields aglucones that rearrange to toxic isothiocyanates. Specifier proteins such as thiocyanate-forming protein from Thlaspi arvense (TaTFP) are non-heme iron proteins, which capture the aglucone to form alternative products, e.g. nitriles or thiocyanates. To resolve the electronic state of the bound iron cofactor in TaTFP, we applied continuous wave electron paramagnetic resonance (CW EPR) spectroscopy at X-and Q-band frequencies (∼9.4 and ∼34 GHz). We found characteristic features of high spin and low spin states of a d 5 electronic configuration and local rhombic symmetry during catalysis. We monitored the oxidation states of bound iron during conversion of allylglucosinolate by myrosinase and TaTFP in presence and absence of supplemented Fe2+. Without added Fe2+, most high spin features of bound Fe3+ were preserved, while different g'-values of the low spin part indicated slight rearrangements in the coordination sphere and/or structural geometry. We also examined involvement of the redox pair Fe3+/Fe2 in samples with supplemented Fe2+. The absence of any EPR signal related to Fe3+ or Fe2+ using an iron-binding deficient TaTFP variant allowed us to conclude that recorded EPR signals originated from the bound iron cofactor.

DNA methylation is enhanced during Cd hyperaccumulation in Noccaea caerulescens ecotype Ganges.

In this study, we assess the DNA damage occurring in response to cadmium (Cd) in the Cd hyperaccumulator Noccaea caerulescens Ganges (GA) vs the non-accumulator and close-relative species Arabidopsis thaliana. At this purpose, the alkaline comet assay was utilized to evaluate the Cd-induced variations in nucleoids and the methy-sens comet assay, and semiquantitative real-time (qRT)-PCR were also performed to associate nucleus variations to possible DNA modifications. Cadmium induced high DNA damages in nuclei of A. thaliana while only a small increase in DNA migration was observed in N. caerulescens GA. In addition, in N. caerulescens GA, CpG DNA methylation increase upon Cd when compared to control condition, along with an increase in the expression of MET1 gene, coding for the DNA-methyltransferase. N. caerulescens GA does not show any oxidative stress under Cd treatment, while A. thaliana Cd-treated plants showed an upregulation of transcripts of the respiratory burst oxidase, accumulation of reactive oxygen species, and enhanced superoxide dismutase activity. These data suggest that epigenetic modifications occur in the N. caerulescens GA exposed to Cd to preserve genome integrity, contributing to Cd tolerance.

Genetic and environmental drivers of large-scale epigenetic variation in Thlaspi arvense.

Natural plant populations often harbour substantial heritable variation in DNA methylation. However, a thorough understanding of the genetic and environmental drivers of this epigenetic variation requires large-scale and high-resolution data, which currently exist only for a few model species. Here, we studied 207 lines of the annual weed Thlaspi arvense (field pennycress), collected across a large latitudinal gradient in Europe and propagated in a common environment. By screening for variation in DNA sequence and DNA methylation using whole-genome (bisulfite) sequencing, we found significant epigenetic population structure across Europe. Average levels of DNA methylation were strongly context-dependent, with highest DNA methylation in CG context, particularly in transposable elements and in intergenic regions. Residual DNA methylation variation within all contexts was associated with genetic variants, which often co-localized with annotated methylation machinery genes but also with new candidates. Variation in DNA methylation was also significantly associated with climate of origin, with methylation levels being lower in colder regions and in more variable climates. Finally, we used variance decomposition to assess genetic versus environmental associations with differentially methylated regions (DMRs). We found that while genetic variation was generally the strongest predictor of DMRs, the strength of environmental associations increased from CG to CHG and CHH, with climate-of-origin as the strongest predictor in about one third of the CHH DMRs. In summary, our data show that natural epigenetic variation in Thlaspi arvense is significantly associated with both DNA sequence and environment of origin, and that the relative importance of the two factors strongly depends on the sequence context of DNA methylation. T. arvense is an emerging biofuel and winter cover crop; our results may hence be relevant for breeding efforts and agricultural practices in the context of rapidly changing environmental conditions.

Genetic dissection of seed characteristics in field pennycress via genome-wide association mapping studies.

Field pennycress (Thlaspi arvense L.) is a new winter annual cash cover crop with high oil content and seed yield, excellent winter hardiness, early maturation, and resistance to most pests and diseases. It provides living cover on fallow croplands between summer seasons, and in doing so reduces nutrient leaching into water sources, mitigates soil erosion, and suppresses weed growth. The first ever genome-wide association study (GWAS) was conducted on a pennycress diversity panel to identify marker trait associations with important seed size and composition related traits. The entire population was phenotyped in three total environments over 2 yr, and seed area, length, width, thousand grain weight, total oil, and total protein were measured post-harvest with specialized high-throughput imaging and near-infrared spectroscopy. Basic unbiased linear prediction values were calculated for each trait. Seed size traits tended to have higher entry mean reliabilities (0.76-0.79) compared with oil content (0.51) and protein content (0.37). Genotyping-by-sequencing identified 33,606 high quality genome-wide single nucleotide polymorphism (SNPs) that were coupled with phenotypic data to perform GWAS for seed area, length, width, thousand grain weight, total oil, and total protein content. Fifty-nine total marker-trait associations were identified revealing genomic regions controlling each trait. The significant SNPs explained 0.06-0.18% of the total variance for that trait in our population. A list of candidate genes was identified based on their functional annotations and characterization in other species. Our results confirm that GWAS is an efficient strategy to identify significant marker-trait associations that can be incorporated into marker-assisted selection pipelines to accelerate pennycress breeding progress.

Field pennycress is an excellent winter annual oilseed that can serve as a cash cover crop. Genotyping-by-sequencing is an effective strategy to genotype pennycress affordably for high-quality genome-wide single nucleotide polymorphisms. First-ever mapping study in a field pennycress association mapping population was conducted. GWAS identified 59 significant marker-trait associations for important quantitative traits. Seed size traits had a higher reliability estimate compared with seed composition traits.

Natural variation and improved genome annotation of the emerging biofuel crop field pennycress (Thlaspi arvense).

The Brassicaceae family comprises more than 3,700 species with a diversity of phenotypic characteristics, including seed oil content and composition. Recently, the global interest in Thlaspi arvense L. (pennycress) has grown as the seed oil composition makes it a suitable source for biodiesel and aviation fuel production. However, many wild traits of this species need to be domesticated to make pennycress ideal for cultivation. Molecular breeding and engineering efforts require the availability of an accurate genome sequence of the species. Here, we describe pennycress genome annotation improvements, using a combination of long- and short-read transcriptome data obtained from RNA derived from embryos of 22 accessions, in addition to public genome and gene expression information. Our analysis identified 27,213 protein-coding genes, as well as on average 6,188 biallelic SNPs. In addition, we used the identified SNPs to evaluate the population structure of our accessions. The data from this analysis support that the accession Ames 32872, originally from Armenia, is highly divergent from the other accessions, while the accessions originating from Canada and the United States cluster together. When we evaluated the likely signatures of natural selection from alternative SNPs, we found 7 candidate genes under likely recent positive selection. These genes are enriched with functions related to amino acid metabolism and lipid biosynthesis and highlight possible future targets for crop improvement efforts in pennycress.

Efficient Protoplast Isolation and DNA Transfection for Winter Oilseed Crops, Pennycress (Thlaspi arvense ) and Camelina (Camelina sativa ).

Pennycress (Thlaspi arvense) and camelina (Camelina sativa) are nonfood winter oilseed crops that have the potential to contribute to sustainable biofuel production. However, undesired agronomic traits of pennycress and camelina currently hinder broad cultivation of these plants in the field. Recently, genome editing using the CRISPR-Cas technology has been applied to improve poor agronomic traits such as the weedy phenotype of pennycress and the oxidation susceptible lipid profile of camelina. In these works, the CRISPR reagents were introduced into the plants using the Agrobacterium-mediated floral dipping method. For accelerated domestication and value improvements of these winter oilseed crops, DNA-free genome editing platform and easy evaluation method of the CRISPR-Cas reagents are highly desirable. Cell wall-free protoplasts are great material to expand the use of gene engineering tools. In this chapter, we present a step-by-step guide to the mesophyll protoplast isolation from in vitro culture-grown pennycress and soil-grown camelina. The protocol also includes procedures for DNA transfection and protoplast viability test using fluorescein diacetate. With this protocol, we can isolate an average of 6 × 106 cells from pennycress and 3 × 106 cells from camelina per gram of fresh leaf tissues. Using a 7.3 kb plasmid DNA carrying green and red fluorescent protein marker genes, we can achieve an average transfection rate of 40% validated by flow cytometry for both plants.

Chromosome-level Thlaspi arvense genome provides new tools for translational research and for a newly domesticated cash cover crop of the cooler climates.

Thlaspi arvense (field pennycress) is being domesticated as a winter annual oilseed crop capable of improving ecosystems and intensifying agricultural productivity without increasing land use. It is a selfing diploid with a short life cycle and is amenable to genetic manipulations, making it an accessible field-based model species for genetics and epigenetics. The availability of a high-quality reference genome is vital for understanding pennycress physiology and for clarifying its evolutionary history within the Brassicaceae. Here, we present a chromosome-level genome assembly of var. MN106-Ref with improved gene annotation and use it to investigate gene structure differences between two accessions (MN108 and Spring32-10) that are highly amenable to genetic transformation. We describe non-coding RNAs, pseudogenes and transposable elements, and highlight tissue-specific expression and methylation patterns. Resequencing of forty wild accessions provided insights into genome-wide genetic variation, and QTL regions were identified for a seedling colour phenotype. Altogether, these data will serve as a tool for pennycress improvement in general and for translational research across the Brassicaceae.

Purification, characterization and antioxidant activity of selenium-containing polysaccharides from pennycress (Thlaspi arvense L.).

Selenium polysaccharides have antioxidant, anti-tumor and other activities. Tea selenium polysaccharides and Lentinan selenium polysaccharides, etc. have been studied. Pennycress is a super-rich selenium plant and there are few studies about pennycress selenium-containing polysaccharides. In this study, the following researches were carried out on pennycress selenium polysaccharides: Selenium-containing polysaccharides were extracted from the leaves of pennycress (Thlaspi arvense L.) using an efficient subcritical water extraction (SWE) process. After purification, two fractions (Se-PPS1 and Se-PPS3) were obtained and subjected to structural identification. The results showed that the molecular weights of Se-PPS1 and Se-PPS3 were 4.2 × 104 Da and 4.5 × 104 Da, respectively. Se-PPS1 is mainly comprised of glucose, galactose, and xylose. Se-PPS3 consists of glucuronic acid, rhamnose, galacturonic acid, glucose, galactose, xylose, and arabinose. Fourier transform infrared spectroscopy (FT-IR) analysis showed that the two fractions had absorption spectra typical of selenium esters and spectral characteristics of polysaccharides. The glycosidic linkages were determined by Gas chromatography-mass spectrometry (GC-MS) and Nuclear magnetic resonance (NMR). The two fractions have (1 â†’ 6)-ß-D- Galp and T-α-D- Glcp configurations. Moreover, in vitro antioxidant activity assays showed that Se-PPS1 and Se-PPS3 exhibited effective radical-scavenging abilities, suggesting they have potential applications as natural antioxidants in food and medicine.

Cruciferous Weeds Do Not Act as Major Reservoirs of Inoculum for Black Rot Outbreaks in New York State.

Cruciferous weeds have been shown to harbor diverse Xanthomonas campestris pathovars, including the agronomically damaging black rot of cabbage pathogen, X. campestris pv. campestris. However, the importance of weeds as inoculum sources for X. campestris pv. campestris outbreaks in New York remains unknown. To determine if cruciferous weeds act as primary reservoirs for X. campestris pv. campestris, fields that were rotating between cabbage or had severe black rot outbreaks were chosen for evaluation. Over a consecutive 3-year period, 148 cruciferous and noncruciferous weed samples were collected at 34 unique sites located across five New York counties. Of the 148 weed samples analyzed, 48 X. campestris isolates were identified, with a subset characterized using multilocus sequence analysis. All X. campestris isolates originated from weeds belonging to the Brassicaceae family, with predominant weed hosts being shepherd's purse (Capsella bursa-pastoris), wild mustard (Sinapis arvensis), yellow rocket (Barbarea vulgaris), and pennycress (Thlaspi arvense). Identifying pathogenic X. campestris weed isolates was rare, with only eight isolates causing brown necrotic leaf spots or typical V-shaped lesions on cabbage. There was no evidence of cabbage-infecting weed isolates persisting in an infected field by overwintering in weed hosts; however, similar cabbage and weed X. campestris haplotypes were identified in the same field during an active black rot outbreak. X. campestris weed isolates are genetically diverse both within and between fields, but our findings indicate that X. campestris weed isolates do not appear to act as primary sources of inoculum for B. oleracea fields in New York.

Cadmium uptake by a hyperaccumulator and three Pennisetum grasses with associated rhizosphere effects.

Pennisetum grasses (P. purpureum Schumach. 'Purple', P. alopecuroides (L.) Spreng. 'Liren' and P. alopecuroides (L.) Spreng. 'Changsui'), and a cadmium (Cd) hyperaccumulator (Thlaspi caerulescens J.Presl & C.Presl), were grown in soil with four Cd addition levels of 0, 2, 20 and 200 mg/kg. Toxicity symptoms were not observed although growth of all plants decreased as Cd addition increased. Shoot bioconcentration factor (BCFS), the translocation factor (TF) and shoot accumulation of Cd for most plants first increased and then declined as Cd concentrations increased. In contrast, the root bioconcentration factor (BCFR) for T. caerulescens declined and root Cd accumulation for T. caerulescens and two P. alopecuroides cultivars increased consistently as Cd levels increased. P. purpureum had the largest biomass with shoot Cd accumulation similar to that of T. caerulescens, despite lower foliar Cd concentration. Although shoot Cd concentrations of two P. alopecuroides cultivars were lower than for P. purpureum, root Cd concentrations were greater. P. purpureum had Cd BCFS and TF (> 1) at 2- and 20-mg/kg Cd addition treatments, similar to T. caerulescens. P. alopecuroides cultivars had Cd BCFR (> 1) and TF (< 1) at all Cd levels. Roots did not affect rhizosphere pH. However, concentrations of acid extractable Cd in rhizosphere soil were lower than those of corresponding non-rhizosphere soil at all Cd levels for T. caerulescens and P. purpureum; T. caerulescens and P. purpureum did not affect less bioavailable Cd fractions. Concentrations of acid extractable Cd in the rhizosphere of the P. alopecuroides cultivars were not reduced at any Cd level. Differences in Cd accumulation among the three Pennisetum grasses were mainly attributable to root biomass and Cd TFs rather than rhizosphere Cd mobility.

Overexpression of ZNT1 and NRAMP4 from the Ni Hyperaccumulator Noccaea caerulescens Population Monte Prinzera in Arabidopsis thaliana Perturbs Fe, Mn, and Ni Accumulation.

Metalliferous soils are characterized by a high content of metal compounds that can hamper plant growth. The pseudometallophyte Noccaea caerulescens is able to grow on metalliferous substrates by implementing both tolerance and accumulation of usually toxic metal ions. Expression of particular transmembrane transporter proteins (e.g., members of the ZIP and NRAMP families) leads to metal tolerance and accumulation, and its comparison between hyperaccumulator N. caerulescens with non-accumulator relatives Arabidopsis thaliana and Thlaspi arvense has deepened our knowledge on mechanisms adopted by plants to survive in metalliferous soils. In this work, two transporters, ZNT1 and NRAMP4, expressed in a serpentinic population of N. caerulescens identified on the Monte Prinzera (Italy) are considered, and their expression has been induced in yeast and in A. thaliana. In the latter, single transgenic lines were crossed to test the effect of the combined over-expression of the two transporters. An enhanced iron and manganese translocation towards the shoot was induced by overexpression of NcZNT1. The combined overexpression of NcZNT1 and NcNRAMP4 did perturb the metal accumulation in plants.

Genome structure and evolution in the cruciferous tribe Thlaspideae (Brassicaceae).

Whole-genome duplications (WGDs) and chromosome rearrangements (CRs) play the key role in driving the diversification and evolution of plant lineages. Although the direct link between WGDs and plant diversification is well documented, relatively few studies focus on the evolutionary significance of CRs. The cruciferous tribe Thlaspideae represents an ideal model system to address the role of large-scale chromosome alterations in genome evolution, as most Thlaspideae species share the same diploid chromosome number (2n = 2x = 14). Here we constructed the genome structure in 12 Thlaspideae species, including field pennycress (Thlaspi arvense) and garlic mustard (Alliaria petiolata). We detected and precisely characterized genus- and species-specific CRs, mostly pericentric inversions, as the main genome-diversifying drivers in the tribe. We reconstructed the structure of seven chromosomes of an ancestral Thlaspideae genome, identified evolutionary stable chromosomes versus chromosomes prone to CRs, estimated the rate of CRs, and uncovered an allohexaploid origin of garlic mustard from diploid taxa closely related to A. petiolata and Parlatoria cakiloidea. Furthermore, we performed detailed bioinformatic analysis of the Thlaspideae repeatomes, and identified repetitive elements applicable as unique species- and genus-specific barcodes and chromosome landmarks. This study deepens our general understanding of the evolutionary role of CRs, particularly pericentric inversions, in plant genome diversification, and provides a robust base for follow-up whole-genome sequencing efforts.

Genomic analysis of field pennycress (Thlaspi arvense) provides insights into mechanisms of adaptation to high elevation.

BACKGROUND: Understanding how organisms evolve and adapt to extreme habitats is of crucial importance in evolutionary ecology. Altitude gradients are an important determinant of the distribution pattern and range of organisms due to distinct climate conditions at different altitudes. High-altitude regions often provide extreme environments including low temperature and oxygen concentration, poor soil, and strong levels of ultraviolet radiation, leading to very few plant species being able to populate elevation ranges greater than 4000 m. Field pennycress (Thlaspi arvense) is a valuable oilseed crop and emerging model plant distributed across an elevation range of nearly 4500 m. Here, we generate an improved genome assembly to understand how this species adapts to such different environments. RESULTS: We sequenced and assembled de novo the chromosome-level pennycress genome of 527.3 Mb encoding 31,596 genes. Phylogenomic analyses based on 2495 single-copy genes revealed that pennycress is closely related to Eutrema salsugineum (estimated divergence 14.32-18.58 Mya), and both species form a sister clade to Schrenkiella parvula and genus Brassica. Field pennycress contains the highest percentage (70.19%) of transposable elements in all reported genomes of Brassicaceae, with the retrotransposon proliferation in the Middle Pleistocene being likely responsible for the expansion of genome size. Moreover, our analysis of 40 field pennycress samples in two high- and two low-elevation populations detected 1,256,971 high-quality single nucleotide polymorphisms. Using three complementary selection tests, we detected 130 candidate naturally selected genes in the Qinghai-Tibet Plateau (QTP) populations, some of which are involved in DNA repair and the ubiquitin system and potential candidates involved in high-altitude adaptation. Notably, we detected a single base mutation causing loss-of-function of the FLOWERING LOCUS C protein, responsible for the transition to early flowering in high-elevation populations. CONCLUSIONS: Our results provide a genome-wide perspective of how plants adapt to distinct environmental conditions across extreme elevation differences and the potential for further follow-up research with extensive data from additional populations and species.

Physiological and Biochemical Effects of Cd Stress in Thlaspi Arvense L-A Non-Accumulator of Metals.

The field pennycress (Thlaspi arvense L.) is a sensitive non-accumulator of metals. The main objective of the present work was to evaluate the effects of Cd stress in T. arvense using a physiological approach. Four-week-old plants were exposed to increasing concentrations of Cd (0, 10, 50 and 100 µM Cd(NO3)2) for 14 days. Plants were harvested at days 0, 3, 7 and 14 to assess the following parameters: shoot and root length and fresh weight; osmolality and conductivity in leaves; chlorophyll a and b contents and photosynthetic efficiency in leaves (young and expanded); water content, Cd accumulation and nutrient imbalances in leaves (young and expanded) and roots. Thlaspi arvense presented a reduction of 27% in the length of roots exposed to 100 µM Cd, and in plants exposed to 50 and 100 µM Cd, a reduction in the fresh weight of both shoots (53 and 58%, respectively) and roots (58 and 63%, respectively) was observed. In general, in plants exposed to the highest concentrations of Cd, a decrease in leaves osmolality and membrane permeability and in water content of both roots and leaves were observed, as well as several nutrient imbalances. As for the photosynthetic parameters, young and expanded leaves responded differently to Cd stress. Exposed young leaves presented no alterations in photosynthetic efficiency or chlorophyll contents, while expanded leaves of 100 µM Cd-treated plants showed a reduction of up to 67% in chlorophyll contents and a concomitant reduction of 30% in photosynthetic efficiency. Overall, Cd-induced senescence in the non-accumulator plant T. arvense exposed to high concentrations of Cd, which was measured as a decrease in several physiological and biochemical parameters and nutrient imbalances. These Cd-induced alterations resulted in lower plant growth, which might have further implications on plant performance at the population level.

Technologies enabling rapid crop improvements for sustainable agriculture: example pennycress (Thlaspi arvense L.).

Growing concerns over food insecurity and ecosystems health related to population growth and climate change have challenged scientists to develop new crops, employing revolutionary technologies in combination with traditional methods. In this review, we discuss the domestication of the oilseed-producing cover crop pennycress, which along with the development of other new crops and improvements to farming practices can provide sustainable solutions to address malnutrition and environmental impacts of production agriculture. We highlight some of the new technologies such as bioinformatics-enabled next-generation sequencing and CRISPR genome editing in combination with traditional mutation breeding that has accelerated pennycress development as a new crop and a potential model system. Furthermore, we provide a brief overview of the technologies that can be integrated for improving pennycress and other crops and the status of pennycress development using these technologies.

Comprehensive benchmarking of software for mapping whole genome bisulfite data: from read alignment to DNA methylation analysis.

Whole genome bisulfite sequencing is currently at the forefront of epigenetic analysis, facilitating the nucleotide-level resolution of 5-methylcytosine (5mC) on a genome-wide scale. Specialized software have been developed to accommodate the unique difficulties in aligning such sequencing reads to a given reference, building on the knowledge acquired from model organisms such as human, or Arabidopsis thaliana. As the field of epigenetics expands its purview to non-model plant species, new challenges arise which bring into question the suitability of previously established tools. Herein, nine short-read aligners are evaluated: Bismark, BS-Seeker2, BSMAP, BWA-meth, ERNE-BS5, GEM3, GSNAP, Last and segemehl. Precision-recall of simulated alignments, in comparison to real sequencing data obtained from three natural accessions, reveals on-balance that BWA-meth and BSMAP are able to make the best use of the data during mapping. The influence of difficult-to-map regions, characterized by deviations in sequencing depth over repeat annotations, is evaluated in terms of the mean absolute deviation of the resulting methylation calls in comparison to a realistic methylome. Downstream methylation analysis is responsive to the handling of multi-mapping reads relative to mapping quality (MAPQ), and potentially susceptible to bias arising from the increased sequence complexity of densely methylated reads.

Spatial genetic and epigenetic structure of Thlaspi arvense (field pennycress) in China.

Thlaspi arvense (field pennycress) is widespread in temperate regions of the northern hemisphere. We estimated the genetic and epigenetic structure of eight T. arvense populations (131 individuals) in China using amplified fragment length polymorphism and methylation-sensitive amplified polymorphism molecular-marker techniques. We detected low diversity at both genetic (mean = 0.03; total = 0.07) and epigenetic (mean = 0.04; total = 0.07) levels, while significant genetic (FST = 0.42, P < 0.001) and epigenetic (FST = 0.32, P < 0.001) divergence was found across the distribution range. Using Mantel testing, we found spatial genetic and epigenetic differentiation, consistent with isolation-by-distance models. We also identified a strong correlation between genetic and epigenetic differentiation (r = 0.7438, P < 0.001), suggesting genetic control of the epigenetic variation. Our results indicate that mating system, natural selection and gene flow events jointly structure spatial patterns of genetic and epigenetic variation. Moreover, epigenetic variation may serve as a basis of natural selection and ecological evolution to enable species to adapt to heterogeneous habitats. Our study provides novel clues for the adaptation of T. arvense.

Functional analysis of ß-ketoacyl-CoA synthase from biofuel feedstock Thlaspi arvense reveals differences in the triacylglycerol biosynthetic pathway among Brassicaceae.

KEY MESSAGE: Differences in FAE1 enzyme affinity for the acyl-CoA substrates, as well as the balance between the different pathways involved in their incorporation to triacylglycerol might be determinant of the different composition of the seed oil in Brassicaceae. Brassicaceae present a great heterogeneity of seed oil and fatty acid composition, accumulating Very Long Chain Fatty Acids with industrial applications. However, the molecular determinants of these differences remain elusive. We have studied the ß-ketoacyl-CoA synthase from the high erucic feedstock Thlaspi arvense (Pennycress). Functional characterization of the Pennycress FAE1 enzyme was performed in two Arabidopsis backgrounds; Col-0, with less than 2.5% of erucic acid in its seed oil and the fae1-1 mutant, deficient in FAE1 activity, that did not accumulate erucic acid. Seed-specific expression of the Pennycress FAE1 gene in Col-0 resulted in a 3 to fourfold increase of erucic acid content in the seed oil. This increase was concomitant with a decrease of eicosenoic acid levels without changes in oleic ones. Interestingly, only small changes in eicosenoic and erucic acid levels occurred when the Pennycress FAE1 gene was expressed in the fae1-1 mutant, with high levels of oleic acid available for elongation, suggesting that the Pennycress FAE1 enzyme showed higher affinity for eicosenoic acid substrates, than for oleic ones in Arabidopsis. Erucic acid was incorporated to triacylglycerol in the transgenic lines without significant changes in their levels in the diacylglycerol fraction, suggesting that erucic acid was preferentially incorporated to triacylglycerol via DGAT1. Expression analysis of FAE1, AtDGAT1, AtLPCAT1 and AtPDAT1 genes in the transgenic lines further supported this conclusion. Differences in FAE1 affinity for the oleic and eicosenoic substrates among Brassicaceae, as well as their incorporation to triacylglycerol might explain the differences in composition of their seed oil.